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We calculate the time-like pion-photon transition form factor and the pion electromagnetic form 
factor up to next-to-leading order (NLO) of the strong coupling constant in the leading-twist kr 
factorization formalism. It is found that the NLO corrections to the magnitude (phase) are lower 
than 30% (30°) for the former, and lower than 25% (10°) for the latter at large invariant mass 
squared > 30 GeV'^ of the virtual photons. The increase of the strong phases with is obtained, 
consistent with the tendency indicated by experimental data. This behavior is attributed to the 
inclusion of parton transverse momenta fcr, implying that the fcr factorization is an appropriate 
framework for analyzing complex time-like form factors. Potential extensions of our formalism to 
two-body and three-body hadronic B meson decays are pointed out. 

PACS numbers: 12.38.Bx, 12.38.Cy, 12.39.St 



I. INTRODUCTION 



I The kr factorization formalism [ll-Q has been applied to next-to-leading-order (NLO) analysis of several space-like 
Cl^, form factors, such as the pion-photon transition form factor 0, the pion electromagnetic (EM) form factor and 
ji^' the B ^ IT transition form factors [l^]- The calculations are nontrivial, because partons off-shell by are considered 
^ , in both QCD quark diagrams and effective diagrams for meson wave functions. The gauge invariance of hard kernels, 
f~| ■ derived from the difference of the above two sets of diagrams, needs to be verified. The regularization of the light-cone 
singularity in the effective diagrams generates double logarithms, which should be summed to all orders. It has been 
found that the NLO corrections, after the above treatments, are negligible in the pion transition form factor, but 
j>! \ reach 30% in the latter two cases. In this Letter we shall extend the NLO kr factorization formalism to the time-like 
00 ' pion transition and EM form factors. 

' One of the widely a.dopted theoretical frameworks for two-body hadronic B meson decays is the perturbative QCD 
C , (PQCD) approach [ll| based on the kx factorization. It has been shown that factorizable contributions to these decays 
can be computed in PQCD without the ambiguity from the end-point singularity. These computations indicated that 
' sizable strong phases are produced from penguin annihilation amplitudes, with which the direct CP asymmetry in the 
■ B^ — >■ K^TT^ decays was successfully predicted. It is then a concern whether PQCD predictions for strong phases are 
stable against radiative corrections. The factorizable penguin annihilation amplitudes involve time-like scalar form 
factors. Before completing NLO calculations for two-body hadronic B meson decays, it is possible to acquire an answer 
to the above concern by studying the time-like pion EM form factor. Besides, the PQCD formalism for three-body 
, B meson decays [l^l has demanded the introduction of two- meson wave functions [isj, whose parametrization also 
rS ' involves time-like form factors associated with various currents. If PQCD results for complex time-like form factors 
, are reliable, a theoretical framework for three-body B meson decays can be constructed. 

NLO corrections to time-like form factors are derived easily from those to space-like ones by suitable analytic 
continuation from —Q^ to , with denoting the momentum transfer squared. We shall present the kx factorization 
formulas for the time-like pion transition and EM form factors up to NLO at leading twist. Following the prescription 
proposed in both the renormalization and factorization scales are set to the virtuality of internal particles. 

With this scale choice, it will be demonstrated that the NLO corrections to the time-like pion transition and EM 
form factors are under control at leading twist. It implies that PQCD predictions for strong phases of factorizable 
annihilation amplitudes in two-body hadronic B meson decays may be stable against radiative corrections. Moreover, 
we observe the increase of the strong phases of the above form factors with Q^, consistent with the tendency indicated 
by experimental data. It will be explained that this behavior is attributed to the inclusion of the parton transverse 
momenta kx- This consistency supports the kx factorization as a potential framework for studying complex time- like 
form factors and three-body B meson decays. 
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FIG. 1: LO quark diagrams for time-like and space-like pion-plioton transition form factors with ® representing the virtual 
photon vertex. The virtuality of the internal quark is labeled explicitly. 



II. PION-PHOTON TRANSITION FORM FACTOR 



In this section we present the leading-twist NLO factorization formula for the time-like pion-photon transition form 
factor. The leading-order (LO) QCD quark diagram describing 7*((7) T^iPi) l{P2) is displayed in Fig. [TJa), where 
the momentum Pi of the pion and the momentum P2 of the outgoing on-shell photon are chosen as 

Pi = (P+,0,Ot), P2 = (0,P2":Ot), Pt=P2=Q/V2, (1) 

with = = {Pi + P2)^ > being the invariant mass squared of the virtual photon 7*. Figure [IJa) leads to the 
LO hard kernel 



(2) 



where A''c = 3 is the number of colors, e is the polarization vector of the outgoing photon, k = (xP{'',0,kT) is the 
momentum carried by the valence quark, If i/\/2N2 is the leading-twist spin projector of the pion, and the subscript 
fi associated with the virtual photon vertex is implicit on the left-hand side. In the previous works on the space-like 
transition form factor 0, [1343' ^^e internal quark remains ofF-shell by (P2 — fc)^ = —{xQ^ + < as indicated in 
Fig.[ljb). For the time- like case, the internal quark may go on mass shell, and an imaginary part is generated in the 
hard kernel according to the principal- value prescription 

= P^T2-^—^ + i-^S{kl-xQ^). (3) 



Fourier transforming Eq. ^ into the impact-parameter b space, we derive the LO pion transition form factor 

F!if\Q^)^ni^^^[dx( bdbMx)exp[-S{x,b,Q,fi)]H^a'HV^Qb), (4) 
o Jo Jo 

with the pion decay constant Z^, the renormalization and factorization scale ^, the Hankel function of the first kind 
Hq^\ and the twist-2 pion distribution amplitude (DA) 0^. Here we shall not consider the potential intrinsic kx 
dependence of the pion wave function [l7| . because its inclusion would introduce additional model dependence, which 
is not the focus of this work. For example, the intrinsic kx dependence has been parameterized into the different 
Gaussian and power forms in . The Sudakov factor sums the double logarithm as\n^ kx to all orders, and 
takes the same expression for both the space-like and time- like form factors [l9[, since it is part of the universal 
meson wave function. For its explicit expression, refer to d, [23,[2l[. Note that Eq. Q can be obtained from the LO 
space-like pion transition form factor in [l^ by substituting mPllp /2 for the Bessel function i^Tp, as a consequence of 
the analytic continuation — — — > [Q^ + ie) in the hard kernel. 

As stated in the Introduction, the NLO hard kernel is derived by taking the difference of the 0{as) quark diagrams 
and the 0{as) effective diagrams for meson wave functions. The ultraviolet divergences in loops are absorbed into 
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the renormalizcd strong coupling constant as (/i) , and the infrared divergences are subtracted by the nonperturbative 
meson wave functions. The above derivation has been demonstrated exphcitly in Q for the space-hke pion transition 
form factor. We repeat a similar calculation for the time-like pion transition factor, and derive the NLO hard kernel^ 



(5) 



with the NLO correction function 

as{iJ.)CF 



47r 



o 1 1 2 I^T 



21nx+ Utt^ 



2[l-iTT- iire {kl - xQ'^)] In 



(6) 



Cp being the color factor. The imaginary pieces proportional to the step function Q are generated from the 0{as) 
quark diagrams. For the evaluation of the 0{as) effective diagrams, we have chosen the direction of the Wilson lines 
the same as in Q in order to respect the universality of the meson wave function. Equation ([5]) can also be achieved 
by substituting {Q^ + ie) for the virtuality of the external photon, and {xQ"^ — k'^ + ie) for the internal quark in 0, 
and then employing the relations ln(— — ie) = In — in and \n{—xQ^ + k^ — ie) ~ In \ xQ^ — fc^ | — i7r&{xQ^ ^ ^t)- 
Fourier transforming Eq. ([5|) to the b space, we arrive at the NLO fc-r factorization formula for the time-like pion 
transition factor 



dx 



6 JO 
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with 



Kj{x,Q ,kT,^l) = -3\n-^ - - In 
17 n^ 
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being the Euler constant. The function 
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where a denotes the order parameter of the Hankel function, comes from the Fourier transformation of ln^(— xQ^+fc^ — 

ie) in Eq. ([6|). For a small argument p = \fxQb 0, its ma gnitude behaves as \uf^'\p)\ - (1/3) In^ p which 
represents a double-logarithmic correction essentially. The perturbative expansion could be improved by summing the 
double logarithm as ln^[x/((5^6^)] in Eq. ([5]), which arises from the Fourier transformation of the term as ln^(|A;|, — 
xQ'^MQ^) in Eq. ([B]). Strictly speaking, it differs from the threshold resummation of as In^ x performed in [23, and 
deserves a separate study. Besides, there is no end-point enhancement involved in the present calculation, so we shall 
not perform the resummation here for simplicity. Equations ^ and ([7]) will be investigated numerically in Sec. IIVI 



III. PION ELECTROMAGNETIC FORM FACTOR 



We then derive the NLO, i.e., 0{p?^) contribution to the time-like pion EM form factor at leading twist. An LO 
quark diagram for the corresponding scattering ^*{_q) — n^(P\) n~{Pi) is depicted in Fig.l^Ja). We choose light-cone 



coordinates, such that the momenta P\ and P2 are parameterized the same as in Eq. ([TJ with 



The valence quark carries the momentum ki 
The LO hard kernel reads 

r(LO) 



{xiPi , 0, kir) and the valence anti-quark carries k2 



'{xi,kiT,X2,k2T,Q ) = i^T^asCp 



XI Tr [i/2^i7M^i 



[xiQ^ 



k2 



- ie){xiX2C 



<-2T\ 



- ie) ' 



= (i^l+-P2)'>0. 
{0,X2P2,'i^2T)- 



(10) 



^ Compared to Q, three effective diagrams for the self-energy corrections to the Wilson lines have been included in Eq. JSjl. 
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FIG. 2: LO quark diagrams for time-like and space-like pion electromagnetic form factors. 



where the denominators {xiQ^ — k^^^) and {xiX2Q^ — \ii-iT + k2Tp) are the virtuality of the internal quark and 
gluon, respectively. The subscript II denotes that the fcr-dependent terms in both the internal quark and giuon 
propagators are retained. When one of the internal particle propagators goes on mass shell, an imaginary part is 
produced according to the principle- value prescription in Eq. (jS]). 

Fourier transforming Eq. (jlOp from the transverse- momentum space (kiT, ^2t) to the impact-parameter space 
(bi, b2), we obtain a double-6 convolution for the LO time- like pion EM form factor (22j 

-^EM^^(Q^) = OAT ^ / dxidx2 / dbidb2bih2as{^i)xi(l).^{xi)(j)^{x2)c^p[~Sii{xi,bi,X2,b2,Q,^j)] 

^^^c Jo Jo 

X H^^\^^Qb2) [Hl:^\^iQbi)Jo{^iQb2)Q{bi - 62) + i?i'^V^Q&2) Jo(V^Qfoi)e(62 - ^i)] , (11) 

with the Bessel function of the first kind Jo, and the Sudakov exponent S'ii(a;i, fei, 2:2, 62, Q, /i) = S{xi,bi,Q, fi) + 
S{x2, &2i Qy fJ-)- The above expression can also be obtained via analytical continuation of the space-like form factor in 
Fig. IHb) to the time-like region. 

The NLO hard kernel for the space-like pion EM form factor has been computed as the difference between the 
one-loop QCD quark diagrams and effective diagrams in Q. To simplify the calculation, the hierarchy xiQ'^ , X2Q'^ 3> 
xiX2Q^ has been postulated, since the fcr factorization applies to processes dominated by small-x contributions. 
Ignoring the transverse momenta of the internal quarks, the LO hard kernel in Eq. pop reduces to 

„(LO), , , r< Tr [j/'2j/'l7M-^l] 

(y"'(a;iX2Q^ - |kiT + k2Tr + le) 
The Fourier transformation of the above expression leads to a single-6 convolution (23j 

Fj;^'^\Q^)^i — -^—^ dxidx2 dbb as{^i)<j)^ {xi)(j)Tr{x2) exp[-5'i(a;i, 3:2, 6, Q, ^i)\H''^\y' XiX2Qb) , (13) 

with the simplified Sudakov exponent Si{xi,X2,b,Q, fi) = Su{xi,b,X2,b,Q, ji). Comparing the outcomes from 
Eqs. (fTTj) and ifTS)) . we can justify the proposed hierarchical relation, and tell which particle propagator, the internal 
quark or the internal gluon, provides the major source of the strong phase. 

Substituting (Q^ + ie) for the virtuality of the external photon, and {xiX2Q^ — |kiT + k2Tp + «e) for the internal 
gluon in 0, we have the NLO hard kernel for the time-like pion EM form factor 

HEM'^\xi,kiT,X2, k2T, Q^, m) = h^M{xi,X2, 5i2,Q, ^i)H'^'^\xi, kiT , X2, k2T, Q^), (14) 
with the NLO correction function 



IEm{Xi,X2,0i2,Q,H) = 



3 u2 17 2Y 13 31 

— m —7 m x-i ~\ mx-i mx2 In Xi H In Xn 
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(15) 
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and the notation 



In Si 



In 



|kiT+k2T| ~XiX2Q'^ 



6(|kiT + k2T| - X1X2Q' 



(16) 



Fourier transforming Eq. (jl4p . we derive the kx factorization formula for the NLO contribution at leading twist 
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with the function 
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(18) 



The perturbative expansion could be im prov ed by organizing the double logarithm as In xi in Eq. (|15p into the 
threshold resummation factor St{xi,Q^) |20l|. This double logarithm, the same as analyzed in j20| . appears in the 
loop correction to the virtual photon vertex under the hierarchical relation xiQ^ ^ Because there is no 

end-point enhancement involved at leading twist, we shall not perform the threshold resummation here. However, 
the end-point enhancement exists in the two-parton twist-3 contribution, for which St will play a crucial role, and be 
implemented in Sec. IIVI We shall investigate the NLO effect at leading twist in the time-like pion EM form factor 
based on Eqs. (HH and p7)) . 



IV. NUMERICAL ANALYSIS 



The numerical analysis is performed in this section, for which we adopt the standard two-loop QCD running coupling 
constant as{^) with the QCD scale Aqcd = 0.2 GeV, the pion decay constant = 0.131 GeV, the nonasymptotic 
two-parton twist-2 pion DA 

(j)^{x) = 6x(l - a;) [1 -I- a2Cf ^(1 - 2x) \ , (19) 

with the Gegenbauer coefficient 02 = 0.2 being fixed by lattice QCD and the Gegenbauer polynomial C|^^(w) = 
(3/2)(5m2- 1). 

We compute the LO and NLO contributions to the time-like pion-photon transition form factor at leading twist via 
Eqs. Q and ([7]), with the renormalization and factorization scale /i being set to the virtuality of the internal quark 
^ = max(Y^(5, 1/&)- The behavior of Q^Ft^^[Q'^) for < 20 GeV^ displayed in Fig. [3] reflects the oscillatory nature 
of the LO hard kernel in the b space. The LO time-like pion transition form factor exhibits an asymptotic magnitude, 
Q2 |i^^.^(Q^)| « 0.225 GeV at large Q^. RecaU that an asymptotic scaling is known as Q'^Ft^^{Q'^) v^/tt = 0.185 
GeV for the space- like pion transition form factor at large —q^ = Q^ [Hj. The larger asymptotic value for the former 
is expected in the kr factorization, because the internal quark may go on mass shell for a time-like momentum transfer 
q, but it is always off-shell for a space-like q. The ratio between the asymptotic values of the above two transition 
form factors is roughly 1.22, comparable to the data 1.14 from the 77 — 7 transition form factors for > 100 GeV^ 
[26| . Note that the time-like and space- like transition form factors would have equal magnitudes in the LO collinear 
factorization without including the parton transverse momentum fc^. The NLO contribution to the time- like pion 
transition form factor is also displayed in Fig. [31 which decreases with as expected in PQCD. Compared to the LO 
result, the NLO correction to the magnitude is about 30% at = 30 GeV^, and less than 20% for > 50 GeV^. 

For the phase, the LO result arises with Q^, and approaches an asymptotic value close to 180° as shown in Fig.|3l It 
is obvious that the variation with is also attributed to the inclusion of the parton transverse momentum kr- If kr 
in Eq.([5]) was dropped, the LO hard kernel reduces to the traditional expression in the collinear factorization, which 
always leads to a real F^^j. A quantitative understanding can be attained via Eq. ([3]): the contributions from the two 
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FIG. 3: Magnitude and phase of the time-hke pion-photon transition form factor at LO (dashed) and up to NLO (solid). The 
NLO correction is marked in gray. 



terms in Eq. Q are comparable at low Q^, such that the time- like pion transition form factor acquires a nontrivial 
phase. At high > 20 GeV^, the phase is dominated by the first term in Eq. ([3]), since it is unlikely to have a large 
parton = xQ^ demanded by the second term. That is, the tiny deviation (less than 5°) of the asymptotic phase 
from 180° is caused by the power-suppressed k^/Q^ effect. The NLO correction to the phase is about 30° at = 30 
GeV^, and fewer than 20° for > 50 GeV^. The above investigation implies that higher-order corrections to the 
complex time-like transition form factors are under control in the kr factorization. As stated before, the perturbative 
expansion could be improved by resumming the double logarithm In^ x in Eq. ([8]) . 

For the analysis of the time-like pion EM form factor, we first identify the major source of the strong phase by 
comparing the results from Eqs. pT|) and in Fig. The renormalization and factorization scale ^ is set to 
/i = max (^y/xiQ, 1/bi, I/62) [S associated with the virtuality of the internal particles. The curve from Eq. PT|) 
implies that the magnitude of the time-like pion EM form factor has an asymptotic behavior Q'^\Fem{Q'^)\ 0.14 
GcV^ as Q 00. Similar to the case of the transition form factor, this asymptotic value is larger than that of the 
space- like pion EM form factor Q, because of the inclusion of the parton transverse momenta kx- The inclusion of 
kx also leads to the variation of the phase with Q^, which arises from the first quadrant, and then approaches an 
asymptotic value close to 165°. For > 15 GeV^, the difference between using the single-6 and double-6 convolutions 
is insignificant in both magnitude and phase, verifying the hierarchical relation xiQ^,X2Q^ 3> xiX2Q'^,k^, and the 
major source for the strong phase as the internal gluon propagator. We then investigate the NLO effect in the time-like 
pion EM form factor based on Eqs. (|lip and ([T7|) . which is also shown in Fig. 2] For > 10 GeV^, the observed NLO 
correction is roughly 25% for the magnitude, and less then 10° for the phase. That is, the perturbative evaluation of 
the time-like pion EM form factor is stable against radiative corrections at leading twist. 



At last, we include another piece of subleading effects, the LO two-parton twist-3 contribution [27[, for completeness. 
Following the same derivation of the twist-2 contribution in Sec. IIIIl the kx factorization formula for the LO two- 
parton twist-3 contribution to the time- like pion EM form factor was obtained in [l^l, where an explicit double-6 
convolution expression similar to Eq. (jll[) can be found. The hard kernel in the impact-parameter space is identical 
to the one in Eq. (fTTj) . We employ the asymptotic two-parton twist-3 DAs, 

^^ix) = l, (^lix)^l-2x, (20) 

with the associated chiral scale fi^^ = 1-3 GeV. The threshold resummation factor St{x,Q) with a shape parameter 
c = 0.4 is included, since the important double logarithm as In^ x at small x needs to be summed (20j . The numerical 
outcomes for the time-like pion EM form factor are presented in Fig.[5l where the available experimental data [28l.[29| 
are displayed for comparison. It is known that the pion EM form factor is dominated by the two-parton twist-3 
contribution, instead of by the twist-2 one at currently accessible energies, because of the end-point enhancement 
developed by the above DAs [2^[33,[3l|. This enhancement is understood easily as follows: the virtual quark and 
gluon propagators behave like 1/xi and l/(a;ia;2), respectively, as indicated in Eq. ([TU]). The twist-2 pion DA is 
proportional to (p-^ix) ^ 0{x), but the twist-3 pion DAs remain constant (^^''^{x) ~ 0(1) at small a;, which then 
enhance the end-point contribution dramatically. This enhancement was also observed in perturbative evaluation of 
the i? — )- TT transition form factors [s^, and confirmed by the light-cone sum- rule analysis [1^. The relative phase 
between the twist-2 and two-parton twist-3 pieces is about 70° as indicated by Figs. [Hand [SJ so the magnitude of the 
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FIG. 4: Magnitude and phase of the time-hke pion EM form factor at leading twist. Contributions from LO with the single-6 
convolution (dotted), LO (dashed), and LO+NLO (solid) are shown. 

form factor is hardly affected by the former. However, tlie twist-2 contribution does have a sizable effect on the phase 
as illustrated in Fig. [S] 

The predictions for the magnitude of the time-like pion EM form factor from the fc^ factorization can accommodate 
the data for > 4 GeV^, an observation consistent with that from the LO analysis [2^. We point out that 
the measured magnitude of the time-like pion EM form factor is larger than the space- like one [22] , and simultaneous 
accommodation of both data is possible in the fcy factorization, but not in the coUinear factorization. Though the 
perturbative calculations may not be justified for small < 4 GeV^, it is interesting to see the coincidence between 
the increases of the phase with from the kr factorization and from the data for < 1.3 GeV^. In a Breit-Wigner 
picture, the observed phase increase could be attributed to a resonant p meson propagator (29j . It happens that the 
parton transverse momentum fc^ plays the role of the p meson mass, such that the two curves in Fig. [5] exhibit the 
similar tendency, and begin to merge for > 1 GeV^. Again, this coincidence cannot be achieved in the collinear 
factorization, which does not generate a significant phase shift. 

The consistency between the present analysis and the data supports the kr factorization formalism as an appropriate 
framework for studying complex time-like form factors. It has been understood that the complex penguin annihilation 
contribution is essential for explaining direct CP asymmetries in two-body hadronic B meson decays This 
contribution involves time-like scalar form factors, which can be calculated in the same fcr factorization formalism. 
It has been observed that the phase of the 5'-wave component in tttt scattering shows a similar dependence to 
that of the P-wave [1^. Therefore, the PQCD predictions for the above direct CP asymmetries are expected to be 
reliable. The formalism for three-body hadronic B meson decays [l^ has required the introduction of two-meson wave 
functions, whose parametrization also involves time-like form factors of various currents. Stimulated by our work, we 
have the confidence on computing these complex time-like form factors directly in the PQCD approach. 

V. CONCLUSIONS 

In this Letter we have calculated the time-like pion-photon transition and EM form factors up to NLO in the fc^ 
factorization formalism. The corresponding NLO hard kernels were derived by analytically continuing the space-like 
ones to the time-like region of the momentum transfer squared Q^. We have identified the fcy-dependent internal 
gluon propagator as the major source for the strong phase of the time-like pion EM form factor, which increases 
with Q^, and approaches an asymptotic value [34| . The magnitudes of the time-like form factors are larger than 
those of the space-like ones. It has been realized that the above features are attributed to the inclusion of the parton 
transverse momenta, and consistent with the tendency implied by the data. It was observed that the NLO corrections 
in magnitude (phase) change the LO leading-twist results by roughly 30% (30°) for the pion transition form factor, 
and 25% (10°) for the pion EM form factor as > 30 GeV^. The stability against radiative corrections justifies 
the kx factorization formalism for both time-like form factors at leading twist. Therefore, the predictions for strong 
phases of annihilation contributions to two-body hadronic B meson decays in the PQCD approach may be reliable. 
The framework presented here will have other applications, for example, to the construction of the two-meson wave 
functions for three-body B meson decays. 
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FIG. 5: Contributions to the time-like pion EM form factor from two-parton twist-3 LO (dotted), two-parton twist-3 LO plus 
twist-2 LO (dashed), and two-parton twist-3 LO plus twist-2 up to NLO (solid). 
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